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Abstract 

This work is an analytical calculation of the path probability for 
random dynamics of mechanical system described by Langevin equa¬ 
tion with Gaussian noise. The result shows an exponential dependence 
of the probability on the action. In the case of non dissipative limit, 
the action is the usual one in mechanics in accordance with the pre¬ 
vious result of numerical simulation of random motion. In the case 
of dissipative motion, the action in the exponent of the exponential 
probability is just the one proposed in a previous work (Q.A. Wang, 
R. Wang, arXiv:1201.6309), an action defined for the total system 
including the moving system and its environment receiving the dissi¬ 
pated energy. In both cases, the result implies that the most probable 
paths are the paths of least action which, in the limit of vanishing 
randomness, become the regular paths minimizing the action. 
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1 Introduction 


A characteristic of random motion is the existence of many possible paths 
between two given positions. A trivial experiment to show this is to drop a 
small feather in the air from a top point to a fixed point below. The first 
time it goes down along a zigzag path, the second time another path, and so 
on. If each path is defined as a band with a reasonable thickness instead of a 
single geometrical line, it is possible to assign a probability of occurrence to 
each of these paths, or the probability for a particle to take. In principle, this 
probability can be experimentally measured by observing a large number of 
identical feathers moving down between the two points, and estimated by the 
number of feathers observed along this path, divided by the total number of 
feathers arriving at the end point below. In practice, this kind of experiments 
is not that easy because, in order to get reliable calculation of probability, 
one has to count a very large number of identical bodies moving between two 
points under identical conditions. 

In a recent work[l], we have proposed numerical measure of this proba¬ 
bility by simulating the random motions subject to conservative forces. The 
aim of this probability measure is to determine its random variables. The 
random motion was generated by a Langevin equation with Gaussian noise. 
This numerical model is different from the above example of feathers in that 
there is no friction or energy dissipation on the falling bodies. This is an 
ideal motion or an asymptotic model for weakly damped random motion. 
The choice of this model, fully discussed in [1], is motivated by a theoretical 
study of path probability in the framework of a probabilistic Hamiltonian 
mechanics based on a generalization of least action principle of the regular 
Hamiltonian mechanics[3, 4, 5]. According to this theory, if the system is 
perturbed by a Gaussian white noise and statistically remains Hamiltonian 
one without or with very weak dissipative force, the path probability should 
be exponentially decreasing function of the action A = f (K — V)dt 
where K is the kinetic energy, V the potential energy, 7 a characteristic 
parameter of the randomness and the time integral is carried out along the 
considered path. This probability is a classical analog of the Feynman’s fac¬ 
tor of quantum Hamiltonian mechanics[2], the difference being that 7 is 
a real number since the path probability must be real. The output of the nu¬ 
merical experiment has verified this theoretical prediction[l] and confirmed 
the generalization of the regular least action principle under these conditions. 

A question following this work is whether or not the same exponential 
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path probability exists for dissipative random motion. The answer to this 
question is crucial for the description of random motions in nature since 
most of them, if not all, are dissipative or over-damped. In this work, a 
analytical derivation of the path probability is provided hrst for the same 
dynamic model as in [1] in the zero dissipation limit and then for the case with 
important dissipation. The reader will hnd the above mentioned exponential 
path probability exists for both models. In the zero dissipation limit, the 
action is just the usual one as expected, while in the presence of important 
dissipation, the action is a new one suggested in a generalization of least 
action principle to dissipative motion[7]. This new action recovers the usual 
action when the dissipation tends to zero. 


2 The model 


The model can be depicted as the following Langevin equation 



dV{x) 

dx 


+ fd + R 


( 1 ) 


where x is the one dimensional position, t the time, V{x) the potential energy, 
fd the friction force, and R the random forces as a white noise. A discrete 
solution of this equation is 


Xi = Xi-i+ri + f{ti) - ( 2 ) 

which is a superposition of a Gaussian distributed random displacement rj 
(solution of the stochastic equation = R) and a regular trajectory /(tj) 
(the solution of the Newtonian equation + fd) corresponding 

to the least action path. Although Eq.(l) is not necessarily a linear equation 
in X, Eq.(2) is always its solution since the acceleration derivative ^ is a 
linear operator. Figure 1 illustrates some sample paths produced with Eq.(2) 
between two positions a and b, each of them being a tube of hnite thickness 
6 (a band in the x — t representation). 

In the numerical experiment of the work [1], the number Nk of the parti¬ 
cles (or geometrical lines) going through a given sample path /c, i.e., all the se¬ 
quences of positions {xo,xi,X 2 ■ ■ -Xn-i^Xn} satisfying {zk,i—d/2 < Xi < Zk,i + 
5/2} for alH = 1 , 2 ,. .., n, is counted, where {a = Zkfl, Zk,i, Zk, 2 , • • •, Zk,n-i, Zk,n 
b} are the axial line of that tube and n the total number of steps. The prob¬ 
ability that the path k is taken is given hy Pk = Nk/N where N is the total 
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Figure 1; Sample paths created with Eq.(2) for a random motion subject to 
the Harmonic force and the Gaussian noise. 
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number of particles (trajectories) moving from a to b through all the consid¬ 
ered sample paths. The discrete normalization of Pk is -Pfc = 1 where s 
is the total number of paths. 

For each sample path, the instantaneous velocity of the particles of mass 
m at the step i is calculated hy Vki = along the axial line. This 

velocity can be approximately considered as the average velocity of all the 
trajectories passing through the tube if 6 is small with respect to the whole 
scale of the motion. The kinetic energy is given by Kk,i = r The action 

along a path k is given by Ak = YA=Ak^'>^k,i “ 

The Gaussian distribution of r* at the step i is given by 


p(,ri,a) 



( 3 ) 


where, in the case of a free Brownian particle for instance, a = 

V^DAt is the standard deviation, D the diffusion constant, and At = fj—fj_i 
is the time interval of each step of the discrete motion. 


3 Path probability of random motion in the 
non dissipative limit 


This is the case where, in the Eq.(l), the dissipative force fd is very weak 
with respect to the conservative force — or fhe dissipative energy by the 
force fd is small compared to the variation of the total mechanical energy. 

From the dehnition of the Wiener measure[6], the probability Pk of a 
given path k can be expressed as follows: 


Pk 


nr=i Pk,i 

J2l=iUtiPk,i 


( 4 ) 


where pk,i = p{zk,i — 5/2 < Xi < Zk,i + 5/2) for all f = 1,2,... ,n, is the 
Gaussian distribution of all the realization Xi of the random variable Xi that 
are found in the tube of the path k. From Eq.(2), pk^i is given by 


Pk,i — 


\/^( 


expj--^- y-T - -}du 




2a2 


( 5 ) 
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For small 5, we have 




\/ 2 ^( 


TTCr 


2 a2 


This leads to 


1 1 " 

Pfc ~ 7 exp{-— '^[zk,i - Xi-i - f{ti) + 


2=1 


where the normalization constant is given by 


Z = ^exp {-—- f{ti) + 
k=l i=l 


This implies 


InPfc 


2 a 2 


'^[Zk,i - Xi-i - f(ti) + - InZ. 


2 = 1 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 


As Z is constant, the function InP^ depends only on the sum over i. Let it 
be Sk, we have 


Sk 


+ '^{Zk,i - Zk,i-lY + '^{Zk,i-l 
2=1 2 = 1 2=1 
n 

^ ^ 2 —l) 

2=1 

n 

+ 2^(zfc,i-i - Xi_i){f{ti) - f{ti-i)) 

2=1 

n 

-2'^{zk,i - Zk,i-i){f{ti) - f{ti-i)) 

2 = 1 


Xi-l) 


2 


( 10 ) 


We remark that : 

• The first term is constant because it is a sum over the deterministic path 

m. 

• The second term in Eq.(5) is just the time sum (time integral in continuous 


version) of the kinetic energy, i.e., 'EY=ii^k,i “ ^k,i-iY = ^ Yh=i 
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• The third term ~ < n6'^ can be neglected for small 6. 

• The fourth term is such that 

n 

2 y^{zk,i - Zk^i_i){zk,i-i - Xi_i) < 26n max \zk,i - Zk,i-i\ • (11) 

k=l,...,s 

It follows that this term can be neglected for small <5 because max \zk,i — Zk,i-i \ 

k=l,...,s ’ ’ 

2=l,...,n 

is constant. 

• The fifth term is such that 

n 

-Xi-i){f{ti) - /(tj-i)) < 25n max - /(ti-i)| • ( 12 ) 

2 = 1 ,...,n 

2=1 

Consequently, this term can be neglected for small 5 because max |/(fi) — /(fj-i) | 

2 = 1,...,22 

is constant. 

• Finally, the sixth term in Eq.(5) reads 


n—1 


-2 '^{zk,i - Zk,i-i){f{ti) - = 2 ^ Zk,i (/(ti-i - 2f{U) + /(ti+i)) 


2=1 


2=1 


+2Zk,o[f{ti) - /(to)] + ‘2Zk,n[f{tn) - /(tn-l)] 
^-1 / f{ti+l)-fiti) _ /(ti)-/hi-l) 

=2j2^k,ii^tf 


At 


At 


i=l 


At 


+2a[f{ti) - /(to)] + 26[/(tn) - /(tn-l)] 

It follows that for (At) small, we have 


(13) 


n—1 


-2^{zk,i- Zk,i-i){f{ti)-f{ti_i)) ^2^Zk,i{Atf (/(tj)) 
2=1 2=1 


+2a[/(ti) — /(to)] + 26[/(t„) — f{tn-i)] 


At this stage, Sk reads 


Sk = 


2At 


m 


n—1 


+ ^k,ijnf{ti) 


2=1 


2 = 1 


(14) 


+ constant (15) 
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where constant includes all the constant terms independent of k mentioned 
above. Since mf{ti) = = -V'{f{ti)), we get 


S,= 


2At 


m 


n—l 


~ Y1 


2=1 


2 = 1 


(16) 


to a constant. In order to introduce the potential energy V{zk,i) on the path 
k, we consider the following approximation 


V{zk,) ^ V{f{U)) + {zk, - f{U))V'{f{U)), (17) 

or Zk,iV\f{ti)) ^ V{zk,i) -V{f{ti)) + f{ti))V'{f{ti)), which leads to 


S,= 


2At 


m 


n—l 




(18) 


_2=1 2=1 

to a constant including all the terms independent of of the path k in the 
Eq.(17). Let us introduce the usual action of the mechanics given by 


n 

Ak = ^Lk,iAt (19) 

2 = 1 

where the Lagrangian is — V(zk^i), the kinetic energy Kk^i = 

and the potential energy V{zk,i), all along the path k, up to the time ti. Fi¬ 
nally, the probability of the path k is given by 

Ffc = ^exp(-?7Afc) (20) 

with T] = All the terms independent of k vanish after the normaliza¬ 
tion over k with the constant Z = 6xp(—r^A^). As r] is always positive, 

Eq.(20) means that the most probable path must be the paths of least ac¬ 
tion Ak- This probability has been observed in the numerical experiment of 
random motion in the weak dissipative limit in the work [1]. 


4 Path probability of dissipative random mo¬ 
tion 

The usual action loses its role of characteristic variable for the paths in the 
variational calculus when dissipation is present. The search for an alternative 







variable has been since long a headache for physicists in analytical mechanics, 
as discussed briefly in [7]. 

In order to introduce the dissipative force fd and the energy it dissipated 
in a Lagrangian to calculate the action, we consider the work Ed of fd, along 
the path /(f*), from the time f = to = 0 to the time t = ti, i.e., Ed{f{ti)) = 
~ lo^' fd{x)dx. Ed, being the energy lost by the system to the environment, 
must be positive. We have, according to the second fundamental theorem of 
calculus, the following expression fd = = ~Ed{f{ti)). 

Now Eq.(16) can be rewritten with the help of the following equation of 
motion 

mfik) = -V'ifiU) - E'.ifit,)) ( 21 ) 

and the approximation 

V(zts) + E^zts) ~ V(fm + Ej(f(tt)) 

+ ( 21 , - f{ti)) [vym + EMiU))] 

After arrangement, Sk reads, to a constant independent of the index k. 


( 22 ) 


S - 


m 


where the action Ak is given by 

n 

^ Lk^iAt 


(23) 


(24) 


2 = 1 


with the Lagrangian j = j — V{zkd) — Ed{zk^i), the kinetic energy is 
Kk,i = \'^Vki, the conservative potential energy V{zk,i) and the energy dis¬ 
sipated Ed{zk,i), along the path k up to the time ti. Finally, the path prob¬ 
ability of the path k is given by 


Pk = -^ewi-i'lAk) 


(26) 


All the terms independent of k vanish after the normalization over k with 
the constant Z = X]fc=i 

The quantity Ak is still called action for dissipative system because, in 
a variational calculus proposed in [7], its vanishing variation 6Ak = 0 can 
lead to the Newtonian equation of motion with dissipative force. It is then 
a characteristic variable of paths for dissipative motion in the same way 
as the usual action for non dissipative motion of Hamiltonian systems. This 
‘dissipative’ action recovers the usual action for the case where Ed is negligible 
with respect to the energy of the system. 
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5 Concluding remarks 

In this work, the path probability for random dynamics of mechanical system 
has been calculated on the basis of a model described by Langevin equation 
with Gaussian noise. The result shows an exponential dependence of the 
probability on the action which conhrms the measure of the path probabil¬ 
ity by numerical simulation of random motion with the same model in the 
non dissipative limit. In the case of dissipative motion, the action can be 
calculated with a Lagrangian which include the energy lost by the system 
due to the dissipation, which conhrms a generalization of the principle of 
least action of Hamiltonian systems to dissipative systems[7]. According to 
the result of this work, for both non dissipative and dissipative systems, the 
exponential path probability implies that the most probable paths are the 
paths of least action which, in the limit of vanishing randomness, become the 
regular paths of least action. 
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